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We propose a quantum dimer model for the metallic state of the hole-doped cuprates at low hole 
density, p. The Hilbert space is spanned by spinless, neutral, bosonic dimers and spin S = 1/2, 
charge +e fermionic dimers. The model realizes a ‘fractionalized Fermi liquid’ with no symmetry¬ 
breaking and small hole pocket Fermi surfaces enclosing a total area determined by p. Exact 
diagonalization, on lattices of sizes up to 8 x 8, shows anisotropic quasiparticle residue around the 
pocket Fermi surfaces. We discuss the relationship to experiments. 


SIGNIFICANCE 

The most interesting states of the copper ox¬ 
ide compounds are not the superconductors with 
high critical temperatures. Instead, the novelty 
lies primarily in the higher temperature metal¬ 
lic “normal” states from which the superconduc¬ 
tors descend. Here we develop a simple, intuitive 
model for the physics of the metal at low carrier 
density, in the “pseudogap” regime. This model 
describes a novel metal which is similar in many 
respects to simple metals like silver. However the 
simple metallic character co-exists with ‘topologi¬ 
cal order’ and long-range quantum entanglement 
previously observed only in exotic insulators or 
fractional quantum Hall states in very high mag¬ 
netic fields. Our model is compatible with many 
recent observations, and we discuss more defini¬ 
tive experimental tests. 

INTRODUCTION 

The recent experimental progress in determining the 
phase diagram of the hole-doped Cu-based high temper¬ 
ature superconductors has highlighted the unusual and 
remarkable properties of the pseudogap (PG) metal — 
see Fig. 1. A characterizing feature of this phase is a 
depletion of the electronic density of states at the Fermi 
energy [1, 2], anisotropically distributed in momentum 
space, that persists up to room temperature. 

Attempts have been made to explain the pseudogap 
metal using thermally fluctuating order parameters; we 
argue below that such approaches are difficult to rec¬ 
oncile with recent transport experiments. Instead, we 
introduce a new microscopic model that realizes an al¬ 
ternative perspective [3], in which the pseudogap metal 
is a finite temperature (T) realization of a novel quan¬ 
tum state: the fractionalized Fermi liquid (FL*). We 



FIG. 1. Schematic phase diagram of hole-doped cuprates 
(apart from those with La doping) as a function of temper¬ 
ature (T) and hole density (p). The antiferromagnetic (AF) 
insulator is present near p = 0, and the d- wave superconduc¬ 
tor (dSC) is present below a critical temperature T c . The 
pseudogap (PG) is present for T < T* and acquires density- 
wave (DW) order at low T. The metallic states are the PG 
metal, the conventional Fermi liquid (FL) and the strange 
metal (SM). The dimer model of the present paper describes 
only the PG metal as a fractionalized Fermi liquid (FL*). 


show that our model is consistent with key features of the 
pseudogap metal observed by both transport and spec¬ 
troscopic probes. 

The crucial observation that motivates our work is 
the tension between photoemission and transport exper¬ 
iments. In the cuprates, the hole density p is convention¬ 
ally measured relative to that of the insulating antiferro- 
magnet (AF), which has one electron per site in the Cu 
d band. Therefore, the hole density relative to a filled 
Cu band, with two electrons per site, is actually 1 + p. 
In fact, photoemission experiments at large hole-doping 
observe a Fermi surface enclosing an area determined by 
the hole density 1 + p [4] , in agreement with the Lut- 
tinger relation. In contrast, in the pseudogap metal a 
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mysterious ‘Fermi arc’ spectrum is observed [5-7], with 
no clear evidence of closed Fermi surfaces. However, de¬ 
spite this unusual spectroscopic feature, transport mea¬ 
surements report vanilla Fermi-liquid properties, but as¬ 
sociated with carrier density p, rather than 1 + p. The 
carrier density of p was indicated directly in Hall mea¬ 
surements [8], while other early experiments indicated 
suppression of the Drude weight [9-11]. While the lat¬ 
ter could be compatible with a carrier density of 1 + p 
but with a suppressed kinetic term, the Hall measure¬ 
ments indicate the suppression of the Drude weight is 
more likely due to a small carrier density. Two recent 
experiments displaying Fermi liquid behavior at low p 
are especially notable: 

• The quasiparticle lifetime r(w,T) determined from 
optical conductivity experiments [12] has the 
Fermi-liquid-like dependence 1/r oc ( hui ) 2 + 
(cnkBT) 2 : with c an order unity constant. 

• The in-plane magnetoresistance of the pseudogap 
[13] is proportional to r _1 (l -f bH 2 r 2 + . ..) in an 
applied field H 1 where r ~ T~ 2 and b is a T- 
independent constant; this is Kohler’s rule for a 
Fermi liquid. 

It is difficult to account for the nearly perfect Fermi- 
liquid-like T dependence in transport properties of the 
pseudogap in a theory in which a large Fermi surface of 
size 1 + p [14] is disrupted by a thermally fluctuating 
order. In such a theory, we expect that transport should 
instead reflect the T dependence of the correlation length 
of the order. 

Moreover, a reasonable candidate for the fluctuating 
order has not yet been identified. The density wave 
(DW) order present at lower temperature in the pseu¬ 
dogap regime has been identified to have a d-form fac¬ 
tor [15-19], and its temperature dependence [20-26] indi¬ 
cates that it is unlikely to be the origin of the pseudogap 
present at higher temperature. Similar considerations 
apply to other fluctuating order models [27] based on AF 
or dSC. 

We are therefore led to an alternative perspective [3], 
in which the pseudogap metal represents a new quan¬ 
tum state which could be stable down to very low T, at 
least for model Hamiltonians not too different from re¬ 
alistic cuprate models. The observed low T DW order 
is then presumed to be an instability of the pseudogap 
metal [28-32] . An early discussion [33] of the pseudogap 
metal proposed a state which was a doped spin liquid 
with ‘spinon’ and ‘holon’ excitations fractionalizing the 
spin and charge of an electron: the spinon carries spin 
S = 1/2 and is charge neutral, while the holon is spinless 
and carries charge +e. However this state is incompati¬ 
ble with the sharp ‘Fermi arc’ photoemission spectrum [7] 
around the diagonals of the Brillouin zone: the spin liq¬ 
uid has no sharp excitations with the quantum number of 


an electron, and so will only produce broad multi-particle 
continua in photoemission. 

Instead, we need a quantum state which has long-lived 
electron-like quasiparticles around a Fermi surface of size 
p, even though such a Fermi surface would violate the 
Luttinger relation of a Fermi liquid. The fractionalized 
Fermi liquid (FL*) [34] fulfills these requirements. 

FRACTIONALIZED FERMI LIQUIDS 

The key to understanding the FL* state is the topolog¬ 
ical nature of the Luttinger relation for the area enclosed 
by the Fermi surface. For the case of a conventional FL 
state, Oshikawa [35] provided a non-perturbative proof of 
the Luttinger relation by placing the system on a torus, 
and computing the response to a single flux quantum 
threaded through one of the holes of the torus. His pri¬ 
mary assumption about the many-body state was that its 
only low energy excitations were fermionic quasiparticles 
around a Fermi surface. This assumption then points to a 
route to obtaining a Fermi surface of a different size [36] : 
we need a metal which, in addition to the quasiparticle 
excitations around the Fermi surface, has global topolog¬ 
ical excitations nearly degenerate with the ground state, 
similar to those found in insulating spin liquids [37, 38]. 
In the context of the doped spin liquids noted earlier, we 
obtain a FL* state when the holon and spinon bind to 
form a fermionic state with spin S = 1/2 and charge +e 
(a possible origin of the binding is the attraction arising 
from the nearest-neighbor hopping), and there is a Fermi 
surface with quasiparticle excitations of this bound state 
[39-41] (other possibilities for the fate of this bound have 
also been discussed [42] ). Such a Fermi surface has long- 
lived electron-like quasiparticles and encloses an area de¬ 
termined by density p, and not 1 +p [41, 43-45], just as 
required by observations in the pseudogap metal. Alter¬ 
natively, a FL* phase can also be obtained from Kondo 
lattice models [46, 47], but we shall not use this here. 

Earlier studies have examined a number of phenomeno¬ 
logical and path integral models of FL* theories of the 
pseudogap [40, 41, 43-45] (and in an ansatz for the pseu¬ 
dogap [48]). These models contain emergent gauge field 
excitations, which are need to provide the global topo¬ 
logical states required to violate the Luttinger relation of 
the FL state. But they also include spurious auxiliary 
particle states which are only approximately projected 
out. The gauge field can undergo a crossover to confine¬ 
ment, but the present models do not keep close track of 
lattice scale Berry phases which control the appearance 
of density wave order in the confining state [49]. Here, we 
propose to overcome these difficulties by a new quantum 
dimer model which can realize a metallic state which is 
a FL*. This should open up studies of the photoemis¬ 
sion spectrum, density wave instabilities, and crossovers 
to confinement at low T in the pseudogap metal. 


3 



FIG. 2. (Color online) A typical dimer configuration identi¬ 
fying a state in the Hilbert space. The blue ellipses are the 
bosons D iv which are spinless and neutral. The green rect¬ 
angles are the fermions Fir, a which carry spin S = 1/2 and 
charge +e. The density of the Fi va dimers is p. 

QUANTUM DIMER MODELS 

Quantum dimer models [50-52] have been powerful 
tools in uncovering the physics of spin liquid phases, and 
of their instabilities to conventional confining phases [53— 
55]. Dimer models of doped spin liquids have also been 
studied [50, 56, 57], but all of these involve doping the 
insulating models by monomers which carry charge +e, 
but no spin. Here we introduce a new route to doping, in 
which the dopants are dimers , carrying both charge and 
spin. 

The Hilbert space of our dimer model is spanned by 
the close-packing coverings of the square lattice with two 
species of dimers (see Fig. 2), with an additional two¬ 
fold spin degeneracy of the second species. It can be 
mapped by an appropriate similarity transform [50] to a 
truncation of the Hilbert space of the t-J model. 

The first species of dimers are bosons, Di V , which re¬ 
side on the link connecting the square lattice site i = 
( i x ,i y ) to the site i+fj, where r) = x = (1, 0) or y = (0,1). 
These are the same as the dimers in the Roklrsar-Kivelson 
(RK) model [50], to which our model reduces at zero 
doping. When connecting to the Hilbert space of the 
t-J model, each boson maps to a pair of electrons in a 
spin-singlet state: 

D] v |0> =► T i7) (cl t 4 +j54 +4 + . it ci 4 )[0)/y2, (1) 

where Ci a is the electron annihilation operator on Cu site 
i with spin a and |0) is the empty state with no 

dimers or electrons. The phase factors Y tr) depend upon 
a gauge choice: for the choice made by RK, T iy = 1 and 
T te = (-1 )*». 

The second species of dimers are fermions , with 
a =f4, which carry spin S = 1/2 and charge +e relative 
to the half-filled insulator, and are present with a density 
p. Each fermionic dimer maps to a bound state of a holon 


and a spinon, which we take to reside on a bonding orbital 
between nearest-neighbor sites: 

4a 1°) =* T ^(4a + 4+fl,a)|0)/V^- (2) 

In a three-band model [58, 59], the state |0) can 
be identified with the S = 1/2 state of a hole delocalized 
over a O site and its two Cu neighbors, considered by 
Emery and Reiter [60, 61]. 

Let us stress our assumption that spinon and holon 
bind not because of confinement but because of a short 
range attraction. Therefore, the bound state (2) can 
break up at an energy cost of order the antiferromagnetic 
exchange, and the holon and spinon appear as gapped, 
free excitations which would contribute two-particle con¬ 
tinuum spectra to photoemission or neutron scattering 
spectra. These fractionalized states can be included in 
our dimer model by expanding the Hilbert space to in¬ 
clude monomers, but we will not do so here because we 
focus on the lowest energy sector. As a consequence, 
there is no monomer Fermi surface [43] in the present 
model of the pseudogap metal. 

The states (1) and (2) are precisely those that domi¬ 
nate in the 2-site dynamical mean field theory (DMFT) 
analysis of the Hubbard model by Ferrero et al. [62] : they 
correspond to the S and 1+ states of Ref. 62 respectively, 
which are shown in their Fig. 15 to be the dominant com¬ 
ponents of the ground state wavefunction at small p (see 
also Ref. 63) . The DMFT analysis captures important as¬ 
pects of pseudogap physics, but with a coarse momentum 
resolution of the Brillouin zone. In DMFT, the states 
on the 2-site cluster interact with a self-consistent envi¬ 
ronment in a mean-field way: the equations have so far 
only been solved at moderate temperatures and the na¬ 
ture of the ultimate ground state at low doping remains 
unclear. Our dimer model is a route to going beyond 
DMFT, and to include the non-trivial entanglement be¬ 
tween these states on different pairs of sites in a non¬ 
mean field manner. The local constraints between dif¬ 
ferent pairs of dimers are accounted for, allowing for the 
emergence of gauge degrees of freedom. 

The original RK model can be mapped to a compact 
U(l) lattice gauge theory [51, 53, 54]. In the doped dimer 
models studied earlier, the monomers then carry U(l) 
gauge charges of ±1 on the two sublattices. By the same 
reasoning, we see that the F va fermions carry no net 
gauge charge, but are instead dipoles under the U(l) field. 

We can now describe our realization of the pseudo¬ 
gap metal. We envisage a state where the confinement 
length scale of the compact U(l) gauge field is large, 
and specifically, larger than the spacing between the F va 
fermions. Then the F^a fermions can move coherently in 
the presence of a dipolar coupling to the gauge fluctu¬ 
ations [41], and they will form Fermi surfaces enclosing 
total area p , thus realizing a FL* state. The confinement 
scale becomes large near the solvable RK point in the RK 
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which preserve the dimer constraint and spin rotation 
invariance. Purely fermionic dimer models with similar 
dimer hopping terms have been considered by Pollmann 
et al. [69]. 


RESULTS 
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FIG. 3. (Color online) Terms in the Hamiltonian Hrk + H 1 . 


model [64, 65], near a Higgs transition to a Z 2 spin liq¬ 
uid induced by allowing for diagonal dimers [55, 66-68], 
or more generally near a deconfined critical point [40]. 
Our approach yields a ‘minimal model’ for realizing FL* 
(which can be a stable, deconfined state in the Z 2 spin 
liquid case), and confinement transitions in metals. 

We present results below for the following Hamilto¬ 
nian, illustrated in Fig. 3, acting on the dimer Hilbert 
space described above 


H 

Hrk 


Hi 


Hrk + Hi + H 2 


E \~ J D ix D i+y,x D iy D i+x,y + 


1 term 


+V D ix D Uy,x D ix D i+y,x + 1 term 

E [ _<1 D \x F Uy,xc* F ixa. D i + y,x + 3 termS 


-U D * . F- F D - 

z i-\-x,y iyot ixa. i-\-y,x 


7 terms 


— ^3 D 


fL„f 


D, 


i-\-x-\-y,x iyot i-\-x-\-y ,xa iy 


7 terms 


—to D}. fI.F 


q / / A n . r A i / . 

^ i-\-2y,x iyot i+2y,xot iy 


D- n + 7 terms 


,( 3 ) 


where the undisplayed terms are generated by operations 
of the square lattice point group on the terms above. The 
first term, Hrk, co-incides with the RK model for the 
undoped dimer model at p = 0. Single fermion hopping 
terms are contained in H i, with hoppings ti which are 
expected to be larger than J. A perturbative estimate of 
the dimer hopping amplitudes ti in terms of electron hop¬ 
ping parameters can be found in the SI Appendix. Note 
that all such terms must preserve the dimer close-packing 
constraint on every site, and we have chosen 3 terms with 
short-range hopping; longer-range hopping terms for the 
fermonic dimers are also possible, but expected to de¬ 
cay with distance and are omitted for simplicity. Finally, 
H 2 allows for interactions between the fermionic dimers, 
with terms of the form 


Ho 


E 


(4 S F t 


i-\-y,xot 


_ pt pt p 

ixa. i-\-y,xp/ iyp i-\-x,yot 


( 4 ) 


We now present results for the dispersion and quasipar¬ 
ticle residue of a single fermion described by Hrk + Hi; 
the interaction terms in H 2 play no role here. At a small 
p , the interactions between the fermionic dimers can be 
treated by a dilute gas expansion in p, while the dominant 
contributions to the quasiparticle dispersion and residue 
arise from the interaction between a single fermion and 
the close-packed sea of bosonic dimers. We computed the 
latter effects by exactly diagonalizing the singe fermion 
Hamiltonian on lattice sizes up to 8 x 8 with periodic 
boundary conditions, with the largest matrix of linear 
size 76,861,458. The RK model has two conserved wind¬ 
ing numbers in a torus geometry, and these conservation 
laws also hold for our model: all results presented here 
are for the case of zero winding numbers. We extend 
these results to non-zero fermion density by interpola¬ 
tion in the SI Appendix. 

Our numerical study explored the dispersion of a sin¬ 
gle fermion over a range of values of the hopping param¬ 
eters. We show in Fig. 4 the dispersion e(k) for a single 
F r)a fermion for hopping parameters obtained by a per¬ 
turbative connection on a t-J model appropriate for the 
cuprates at the RK point V = J = 1. The SI Appendix 
has similar results for additional parameter values. 

The minima of the fermion dispersion were found 
at different points in the Brillouin zone, but there 
was a wide regime with minima near momenta k = 
(±7t/2, ±7t/2). In fact, for the momentum points allowed 
on a 8 x 8 lattice, the global minimum of the dispersion 
in Fig. 4 is exactly at (7t/2,7t/2). However, it is also 
clear from the figure that the dispersion is not symmet¬ 
ric about the antiferromagnetic Brillouin zone boundary, 
and that any interpolating function will actually have 
a minimum at {k m ,k m ) with k m < tt/2. A dispersion 
with these properties is of experimental interest because 
it will lead to formation of hole pockets near the min¬ 
ima for the dimer model with a non-zero density of F^ a 
fermions. Fig. 6 shows that changes to the dispersion 
from a 6 x 6 lattice are smaller than 5%. 

Our numerical results also yield interesting informa¬ 
tion on the quasiparticle residue of the electron operator. 
This is non-trivial even for the case of a single fermionic 
dimer, because unlike a free electron, a fermionic dimer 
can only move by ‘resonating’ with the background of 
bosonic dimers, as is clear from Fig. 3. In the presence of 
a finite density of fermionic dimers, there will be an addi¬ 
tional renormalization from the interaction between the 
fermions which we will not compute here. We don’t ex- 
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FIG. 4. (Color online) Lowest energy of a single charge +e 
F va fermion as a function of momentum k. We take hopping 
parameters obtained from the t-J model, t\ = —1.05, £2 = 
1.95, and £3 = —0.6, at the RK point F = J = lona 8 x 8 
lattice with periodic boundary conditions and zero winding 
numbers. Note that the dispersion is not symmetric about 
the magnetic Brillouin zone boundary i.e. across the line 
connecting ( 7 r, 0) to (0, 7 r). Line cuts of this dispersion are in 
the SI Appendix. 



FIG. 5. (Color online) Quasiparticle residue of a charge +e 
fermion computed from Eq. ( 6 ) for the parameters in Fig. 4, 
for a 8 x 8 lattice. The symmetry of the wavefunction yields 
Z{ k) = 0 for all points between (7r, tt) and (-7r,0). Line cuts 
of Z( k) are in the SI Appendix. 


pect this to have a significant k dependence around the 
hole pockets. In the dimer model subspace defined by 
the states in Eqs. (1) and (2), the electron annihilation 
operator on site i has the same matrix elements as 




^a/3 


^ixB^ix 


. F' n „ 

i—x,xBi—x,x 


~^F iy pD iy + F i _. y pD i _~ y S j , ( 5 ) 


relating the site to the 4 bonds around it (e a p is the unit 
antisymmetric tensor). Then the quasiparticle residue is 


e(k) 




FIG. 6 . Line cuts of the dispersion in Fig. 4 (top), and of 
the quasiparticle residue in Fig. 5 (bottom). Also shown are 
the results from exact diagonalization on a 6 x 6 lattice for 
comparison (red squares), which has a different set of allowed 
momentum points. The overall shape of the dispersion re¬ 
mains the same as for the 8 x 8 lattice, and the fractional 
changes to e(k) are smaller than 5%. The inset shows the 
residue between ( 0 , 0 ) and (tv, tv) on a logarithmic scale. 


obtained by computing 

Z(k) = |(^(k)|C a (-k)|T RK >| 2 , ( 6 ) 

where |W RK ) is the ground state of the undoped model 
TIrk, and |\I/jc(k)) is the ground state of L/rk + H\ +H 2 
in the sector with one F^a fermion and total momentum 
k (the energy difference between these two states is e(k)). 
We show the values of Z( k) in Fig. 5, with parameters 
the same as those in Fig. 4. Note the strong suppression 
of the residue in the second antiferromagnetic Brillouin 
zone; line-cut plots of Z( k) in Fig. 6 highlight this sup¬ 
pression. We found this suppression of Z(k) to be a ro¬ 
bust property in the regime of hopping parameters (with 
£2 > 0 ) which had minima in the fermion dispersion along 
the Brillouin zone diagonal. This result implies that the 
quasiparticle residue will be highly anisotropic around 
the hole pockets that appear in the finite fermion density 
case, with little spectral weight along the ‘back side’ of 
the pocket. 

It is also possible to study the system in perturba¬ 
tion theory in ti/J. We begin with the model with one 
fermion at t\ = t 2 = £3 = 0. The problem reduces to 
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that of finding the ground state of Hrk in the presence 
of a stationary fermionic dimer: it is possible to do this 
analytically at the solvable RK point V = <7, as described 
in the SI Appendix. The fermion hoppings at non-zero 
U is then computed perturbatively in a single-particle 
tight-binding model, with hopping matrix elements de¬ 
termined by overlaps of the wavefunctions with a sta¬ 
tionary fermion. At the RK point, each matrix element 
reduces to the evaluation of a dimer correlation function 
in the classical problem of close-packed dimers on the 
square lattice [70]. The computation of these matrix ele¬ 
ments, and the resulting fermion dispersion is described 
in the SI Appendix. This perturbative dispersion is found 
to be in good agreement with our exact diagonalization 
results only for \ti/J\ < 0.01. This rather small upper 
limit is likely a consequence of the gaplessness of the RK 
point, so that higher order corrections involve non-integer 
powers of t t / J. 

DISCUSSION 

In this article we developed a new class of doped dimer 
models featuring coherent electronic quasiparticle excita¬ 
tions on top of a spin-liquid ground state. The scenario 
considered here is based on the assumption that spinons 
and holons form bound states on nearest neighbor sites. 
These fermionic bound states with spin 5 = 1/2 and 
charge +e form a Fermi sea with density p and are ob¬ 
servable as electronic quasiparticles in experiments. Such 
a Fermi sea realizes a topological quantum state called 
the ‘fractionalized Fermi liquid’ [34], whose Fermi sur¬ 
faces encloses an area distinct from the Luttinger value 
in a conventional Fermi liquid. 

The undoped RK model on the square lattice features 
a deconfined spin liquid ground state only at the special 
RK point J = V, whereas the ground state breaks lattice 
symmetries away from this point. Consequently, our nu¬ 
merical results of a single fermionic dimer coupled to the 
background of bosonic dimers focused at the RK point to 
uncover properties of the FL* state. At finite densities of 
fermionic dimers we expect that our model (3) features 
a FL* phase in an extended parameter range. However, 
similar to the RK model, we also expect a wide param¬ 
eter regime where our model has a ground state with 
broken lattice symmetries. We leave the computation of 
the phase diagram of our model for future study. 

The main implication of our model of the pseudogap 
metal (in zero applied magnetic field and at moderate 
T below T*) is that there are 4 well-formed pockets of 
charge +e fermions carrying spin S = 1/2 in the vicinity 
of (but not exactly centered at) momentum ( 7 t / 2 , 7 t / 2 ) . 
The total area enclosed by these pockets is 2tt 2 p. Clearly, 
such pockets can immediately explain the Fermi liquid¬ 
like transport observed in recent optical [12] and magne¬ 
toresistance [13] measurements. We also note that the 


hopping of electrons between CuC >2 layers requires to 
break either fermonic or bosonic dimers in our model, 
which naturally accounts for the observed gap in c-axis 
optical conductivity. 

Experiments which involve removing one electron from 
the system (such as photoemission) have difficulty ob¬ 
serving the ‘back sides’ of the pockets because of the 
small (but non-zero) quasiparticle residue Z(k) noted 
above (see Fig. 5). We propose this feature as an ex¬ 
planation for the photoemission observation of Ref. 7 in 
the pseudogap metal. For further studies of these pock¬ 
ets, it would be useful to employ experimental probes 
of the Fermi surface which keep the electrons within the 
sample [71]: possibilities include ultrasound attenuation, 
optical Hall, and Friedel oscillations. 

Our theory can be loosely summarized by ‘the electron 
becomes a dimer in the pseudogap metal’, as in Eq. (5): 
with a spin-liquid background present, there can be no 
single site state representing an electron, and a dimer is 
the simplest possibility. 

The main advantage of our quantum dimer model 
over previous treatments [40, 41, 43-45] of fractional¬ 
ized Fermi liquids (FL*) is that it properly captures lat¬ 
tice scale dispersions, quasiparticle residues, and Berry 
phases: all of these are expected to play crucial roles 
in the crossovers to confinement and associated symme¬ 
try breaking at low T [49, 53, 54], Given the elongated 
dimer and dipolar nature of the electron, Ising-nematic 
order [72] is a likely possibility; the d-form factor den¬ 
sity wave [15, 16] is then a plausible instability of such 
a nematic metal. The interplay between the monopole- 
induced crossovers to confinement [53, 54] and the density 
wave instabilities of the hole pockets [31, 32] can also be 
examined in such dimer models. The onset of supercon¬ 
ductivity will likely require additional states, such as a 
spinless, charge +2e boson consisting of a pair of empty 
sites. 
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CONNECTION TO THE t-J MODEL 

A quantitive connection between the dimer model Hamiltonian and the t-J model can only be made in the limiting 
regime of a certain large N SU(A') antiferromagnet [52]; we also require t <C J to allow restriction to the dimer 
subspace in the doped case. Eventually, we hope that it will be possible to use e.g. DMFT methods [62] to make this 
connection. Nevertheless, it is useful, for now, to determine the relationship between the hopping parameters of the 
dimer and t-J model in the perturbative approach of large N and small t. We will directly work with the physical 
SU(2) case here. 

The dimer hopping amplitudes are given by the matrix elements 

*1 = -(OltfolO) 

*2 = -<Oltfo|Q> 

h = -{ESItfoltfS . (SI) 

where the fermonic (bosonic) dimer is represented by the black (open) ellipse and Hq denotes the standard tight- 
binding hamiltonian of electrons including up to third-nearest neighbour hopping 

tf0 = -*£ (cLc#a + hx -) - ^2 ( C L C ja + h - c -) - t" (4a C 7 « + h ' C ') ' ( S2 ) 

(ij) ««» «<ii»> 

Using the representation in terms of the dimer operators 

iO) = ^U,4ti°) 

IQ> = ^Wl°> 

IQ> = ^+fi,x4tl°> 

\m= ^i+x+y,x 

iro = Diy ^i+x+y,x f |0) , (S3) 

we can obtain expressions in terms of the electron operators via the mappings in Eqs. (1) and (2). Evaluating the 
matrix elements is now straightforward, and we obtain 

U = ~2 & + *0 

t2 = \(t- t') 

t3 =—-(t + t 1 + t") , (S4) 

where we used the mappings for the dimer operators in Eqs. (1) and (2). Choosing electron hopping parameters 
suitable for cuprates, t'/t = —0.3, t"/t = 0.1, we obtain 


ti = —0.35t 
t2 = 0.65 1 

t 3 = —0.2 1 . (S5) 


The dispersion and residue with the hoppings obeying Eq. (S5) with t = 3 were shown in Figs. 4 and 5 in the main 
text, and their line cuts appear in Fig. 6. 
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FIG. SI. Dispersion and quasiparticles residue for 


e(k) 



= V = 1, and ti = t 2 = —tz = 1 on a 8 x 8 lattice. 


Z( k) 



FIG. S2. Line cuts of the dispersion and quasiparticles residue in Fig. SI, along with corresponding data for a 6 x 6 lattice. 


EXACT DIAGONALIZATION 

This section contains additional plots of the energy eigenvalues and quasiparticle residues obtained by exact diag- 
onalization of H^k + Hi. We note that the exactly solvable RK point of the undoped mode at V = J = 1 can be 
extended to the doped case e.g. by including a nearest-neighbor repulsion V\ = t\ between the F va and D ri dimers, 
and setting t 2 = £3 = 0. However, the exact solution only yields the energy at zero momentum, and numerics are 
required for the full dispersion even at this point. 

In Figs. SI and S2 we show the results for J = V = 1 and t\ = t 2 = —t 3 = 1. We also consider the effect of 
uniformly scaling all the t/ s while keeping J fixed. Figs. S3 and S4 contain results for t\ = t 2 = — 1 3 = 0.5, while 
Figs. S5 and S 6 contain results for ti = t 2 = —t 3 = 2.0. Note that the while the general forms of the dispersion 
and quasiparticle residue remain the same when we uniformly scale the ti, the finite size corrections become larger at 
larger ti/ J. 

Lastly, in Fig. S7, we consider a set of parameters for which the minimum of the dispersion, e(k) remains on the 
diagonal near (7 t/2,7t/ 2), but the quasiparticle residue Z( k) vanishes for k along the diagonal. This happens because 
the wavefunction of the F xa fermions has the opposite sign of the wavefunction for the F ya fermions; then by Eq. (5), 
Z{ k) will vanish for k along the diagonal. 

As an alternative way of presenting the above data, we used the information on the dispersion, e(k), from Fig. 4, and 
the quasiparticle residue, Z( k), from Fig. 5, and interpolated to all k in the Brillouin zone. The resulting electronic 
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- 1.6 

-1.7 

- 1.8 

-1.9 

- 2.0 


FIG. S3. As in Fig. SI, with the parameters ti = t 2 = — = 0.5 at the RK point V = J = 1. 



FIG. S4. As in Fig. S2, with the parameters ti = fe = — £3 = 0.5 at the RK point V = J = 1. 



FIG. S5. As in Fig. SI, with the parameters ti = t 2 = —tz = 2 at the RK point V = J = 1. 
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FIG. S6. As in Fig. S2, with the parameters t\ = = —t^, = 2 at the RK point V = ,7 = 1. 



FIG. S7. Dispersion and quasiparticles residue for J = V = 1, and hopping parameters t\ = —0.5, t 2 = —1 and t 3 = 0.5 on a 
6x6 lattice. For this case, the quasiparticle residue Z( k) vanishes for k along the diagonal because the wavefunction is odd 
under reflections about the diagonal. 


spectral weights at the Fermi energy are shown in Fig. S8 at two different doping densities. 
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FIG. S8. Electronic spectral functions A(k) = Z(k) <5(e(k) — /r) at the Fermi level obtained by interpolating the data in Figs. 4 
and 5. For illustrative purposes, we replace the delta function by a Lorentzian with finite width. Chemical potentials /i are 
chosen to obtain doping densities of p = 0.083 (left) and 0.225 (right). 






6 


PERTURBATIVE COMPUTATION 


This supplement will describe the perturbative evaluation of the single fermion dispersion for small tj under the 
Hamiltonian ILrj< + Hi in Eq. (3) at the RK point V = J. 

At zeroth order in the t,, we need the ground state of a stationary fermion. This ground state is given by an equal 
weight superposition of all dimer configurations with the fermionic dimer in a fixed position: 




(S 6 ) 


where the fermionic dimer is represented by the red bar. The position of the fermionic dimer can be arbitrary, and 
so there is a large degeneracy in the ground state subspace. 

Let us denote the fermionic hopping terms Tf, T 2 , X3 with matrix elements t\, £ 2 , t 3 respectively. These operators 
Ti,T 2 ,T 3 lift the degeneracy of the ground state subspace. Within this subspace, the only non-zero matrix elements 
which are inequivalent under translations and rotations are: 


Zi = 





(S7) 


If we introduce a partition function Q that counts the number of single-flavor dimer close packing configurations 
with a certain number of dimers fixed, we have 


Z\ — Z 2 — 


Q 


Q 



By using a Grassmannian representation of the partition function Q we find 


Z\ — Z 2 — 


1 

2 ’ 



(S 8 ) 


(S9) 


At this point, the problem has been reduced to a “tight-binding model” of a single fermion hopping on the links of 
the square lattice. So the fermion dispersion and wavefunction is determined by diagonalizing a 2 x 2 matrix at each 
momentum k: 


H(k) 


— 2tiZi cos(ky) 

C*(k) 


C(k) \ 
—2tiZicos(k x ) ) 


(S10) 


with 


C(k) = -t 2 Z 2 ( 1 + e ifel )(l + e~ iky ) - t 3 Z 3 [(1 + e lk *){e iky + e~ 2lky ) + (1 + e ~ ik y)(e 2ik * + e~ ik *)\ (Sll) 

We compare this perturbative computation of the quasiparticle energy with the exact diagonalization in Fig. S9 and 
find good agreement at \ti/J\ = 0.01. At larger values of |tj/J| the perturbative computation deviates significantly 
from the numerics, as shown in Fig. S10. 


Details of the Grassmannian approach 


We use Grassmann variables to compute the function Q [70]. The partition function with no restriction can be 
written as 


Q 


= / d?7d?7 e 


f; e s( v ,v) 


(S12) 


where there is a pair (777/) of Grassmann variables for every site of the lattice, and the action is quadratic. It contains 
a pair (77,77) for every colored bond in fig. Sll, with unit coefficient. The circle stands for 77 and the cross for 77, and 
the arrow determines which of the two variables comes first. 

When expanding the exponential, only those strings of 77 s and fjs that contain exactly one pair for every site will 
contribute 1 or —1 to the integral. Pairs marked in red in fig. Sll commute with those marked in blue, and every 
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(7T,7T) k (7T,0) 


( 0 , 0 ): 


FIG. S9. The left panel plots the dispersion obtained from diagonalizing the 2x2 matrix in Eq. (S10), obtained from 
perturbation theory. The right panel compares the perturbative results (shown as the full line) with the exact diagonalization 
of -ffiiK + Hi at t\ — i ‘2 = — f 3 = 0.01 and V = J = 1. The latter results are for lattices of size 6 x 6 (red squares) and 8 x 8 
(blue circles). 



FIG. S10. As in Fig. S9, but at ti = t 2 = — £3 = 0.3 and V = J = 1. The exact diagonalization results here are only for the 
6 x 6 lattice (red squares). 


single red pair must appear in the string. On the other hand, only a subset of the blue pairs can appear in the string, 
and it must form a closed packing dimer covering of the lattice, or else the string will integrate to zero. The non-trivial 
task is to show that the sign of the contribution will be positive. This is achieved by careful choice of the ordering 
structure shown by the arrows in fig. Sll. 
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FIG. Sll. Graphical representation of all terms in the action (S13). The circle stands for q and the cross for fj, and the arrow 
determines which of the two variables comes first. The bonds emphasized by a darker color form a unit cell that covers the 
lattice. 


Introducing explicit integer coordinates for the lattice we have: 



m,n 


^2m,2n+l 

V]2m,2n 

+ V2m+l,2n V2m+l,2n+l + 

(S13) 

^ l2m,2n 

W2m+l,2n 

+ V2m,2n+1 V2m+l,2n+l + 

(S14) 

fj2m+l,2n 

V2m+2,2n 

+ T]2m+l,2n+lV2m+2,2n+l + 

(S15) 

T]2m,2n 

^?2m,2n+l 

+ t?2m+l,2«+lt?2m+l,2n + 

(S16) 

^?2m,2n+l 

^?2m,2n+2 

+ f?2m+l,2rt+2t?2m+l,2n+l • 

(S17) 


In order to simplify the action, we go to momentum space: 

P2m^2n — 

f 4 P 4q e i ( pm+qn ^ipi pq 

V2m+l,2n+l = 

[ 4p4q e^ pm+qn ^ 2pq 

V2m+l,2n 

[ ApM e^ m+qn ^ 3pq 

2 r ■ = 

f dpM e i{pm+qn) ^4 P q 


1]2m,2n 

f 4p4q e~ i(prn+qn) ^ lpq 

(SIS) 

V2m+l,2n+l = 

f &PM e-^ pm+qn ^ 2pq 

(S19) 

V2m+l,2n = 

j 4p4q e~^ pm+qn ^ 3pq 

(S20) 

f?2m,2n+l = 

j 4p4q e- i(pm+qn) ^ ipq 

(S21) 


and we have 


S = J 4p4q(l/j lpq ljj 2 pq i>3pq ^4 pq) 


/0 0 
0 0 
0 1 
\1 0 


1 — e lp 
1 - e~ iq 
0 
0 


1 — e iq 
— (1 — e~ ip ) 
0 
0 


\ 

't 1 


1p2pq 


■03 pq 

/ 

^ipq) 


(S22) 


We may now use the Grassnrannian representation to compute ratios of Q. 
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Normalization 


We know by symmetry that 

1 

4’ 


but it is instructive to obtain this results using the Grassmann approach. Still with reference to fig. Sll 
discussion in the previous section, we have 

Aq — (^?2m,2n^72m,2n-{-l} ? 

where we have chosen the lower left vertical bond, but any other would have worked. We have 

A 0 = jApAqArAs (ipi pq 'tp 4 rs) 

1 - e~ iq _ 1 
2(2 — cosp — cosg) 4 

Here we used the fact that for a quadratic action S = ipKi} we have = [K _1 ]ji■ 



Q 

A 0 = — 

Q 


Z\ and Zo 


We compute 


A\ = 


Q 


Q 


We have 


dll — {V2m,2n r l2m,2n+l r l2m+l,2n+l' r l2m+l,2n) 

= (V2m,2nV2m,2n+l) ( V2m+l,2n+lV2m+l,2n } + (ll2m,2nV2m+l,2n) (V2m,2n+lV2m+l,2n+l) 


11111 

= - X-1-X - = - . 

4 4 4 4 8 


Therefore we have 


Zl = Z2 = Air 2 


Zz 


We compute 


We have 


Aa — 


Q 


Q 


A 3 — {T) 2 m, 2 nV 2 m+l, 2 n'n 2 m+ 2 , 2 n'n 2 m+ 2 , 2 n+l) 

= ( ? 72m,2n^2m+l,2ra) {V 2 m+ 2 , 2 n‘>j 2 m+ 2 , 2 n+l) 


{ r ]2m,2nV2m+2,2n+l) (V2m+l,2nV2m+2,2n) 


1 1 

~ 4 X 4 
_ 1 1 
“ 4 X 4 


ApAq 


e~ ip (l - e~ iq ) 
2(2 — cosp — cos q) 

1 


1 _ 1\ x I = _ 

7T 4 / 4 47T ’ 


(523) 
and to the 

(524) 


(S25) 


(S26) 


(S27) 


(S28) 


(S29) 


(S30) 
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and hence 


^3 


A3 

Ao 


1 

7T 


(S31) 



